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Abstract. We report the results of a search for axions from the 14.4 keV Ml transition from 
57 Fe in the core of the sun using the axio-electric effect in Te02 bolometers. The detectors 
are 5x5x5 cm 3 crystals operated at about 10 mK in a facility used to test bolometers for 
the CUORE experiment at the Laboratori Nazionali del Gran Sasso in Italy. An analysis of 
43.65 kg-d of data was made using a newly developed low energy trigger which was optimized 
to reduce the detectors energy threshold. An upper limit of 0.63 c-kg _1 -d _1 was established 
at 95% C.L.. From this value, a lower bound at 95% C.L. was placed on the Peccei-Quinn 
energy scale of f a > 0.76 x 10 6 GeV for a value of S=0.55 for the flavor-singlet axial vector 
matrix element. Bounds are given for the interval 0.15 < S < 0.55. 
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1 1 Introduction 

2 Quantum chromodynamics or QCD, largely accepted as the best theory describing strong 

3 interactions, contains one curious blemish known as "the strong CP problem". QCD pre- 

4 diets a large neutron electric dipole moment, of the order \d n \ ~ 10~ 16 e-cm, whereas the 

5 experimental bound is |cZ n | < 2.9 x 1CP 26 e-cm [1]. This fact puts an unnaturally small upper 

6 limit (< 10~ 10 ) to the Oqcd parameter, the strength of the CP violating term present in 

7 the QCD. In order to explain this small value Roberto Peccei and Helen Quinn proposed [2] 

8 that the QCD Lagrangian possessed an additional global U(l) symmetry which modified the 

9 CP-violating term to: 

cl 

£e = (Oqcd - y)7^2G% u G afJ ,u- (1.1) 

10 where g s is the strong coupling constant, Ga" the gluon field, a a new scalar field, and f a an 

11 energy scale. Non-perturbative effects induce a potential for the field a that has a minimum 

12 at a = f a 0QCD which causes the spontaneous breaking of the global U(l) symmetry. As 

13 later shown by Weinberg and Wilczek [3], the spontaneous breaking symmetry produces a 

14 Nambu-Goldstone boson known as the axion. 

15 The "standard" Peccei-Quinn axion with symmetry breaking scale of the order of the 

16 electro-weak scale is ruled out by experiments. However, other models of "invisible" axions 

17 which break the symmetry at much higher energies are still an interesting possibility, and 
is the search of an axion mass in the range between 10" 10 - 10" 6 eV is a very active field 

19 of research. The possibility that the axion might be most or part of the dark-matter has 

20 reinforced even further the interest for this field [4]. 

21 The most common models are the so called hadronic or KSVZ axion model [5] and the 

22 GUT or DFSZ models [6]. In both cases, the mass of the axion is directly related to the 

23 axion's decay constant, f a , the pion's mass 771^=135 MeV, the pion's decay constant f n « 92 

24 MeV, and the up-down quark mass ratio, z=0.56 [7], as follows: 

ma= {(l + z + w)(l + z) 

25 where w=m u /m s = 0.029. 
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Figure 1. Axio-electric cross sections of TeC>2 in the energy range 0-20 keV using / a =8.38x 10 6 GeV. 

26 The detection mechanism used in the present work is the axio-electric effect, which is 

27 the equivalent of a photo-electric effect with the absorption of an axion instead of a photon. 

28 The cross section, which can be seen in Fig. 1 for Te02, is given by [8]: 



\m P c 2 J p 



with, a a 



2a E M ym e c 2 ' '" 
1 /2x' e m e c 2 \2 

fa 



(1.3) 
(1.4) 



and a pe is the photo- 



29 For Eq. 1.4, x' e « 1, m e c is the electron mass in GeV, cxem 

30 electric cross sections for each element (taken from [9]). 

31 The axion source in this search is the Ml transition produced by thermal excitation of 

32 57 Fe in the solar core. The isotope 57 Fe is stable and has 2.12% natural abundance, yielding 

33 an average Fe density in the Sun's core of (9.0 db 1.2) 10 19 cm -3 . The uncertainty is mostly 

34 due to the different metal diffusion models in the core and is computed in [10]. The first 

35 excited state is at 14.4 keV, low enough to be thermally excited in the interior of the sun, 

36 which has an average temperature kT ~ 1.3 keV [11] [12]. 

37 In this paper we rely on [12] for the determination of the expected axion flux and axion 

38 interaction rate in Te02 crystals. The 15% error in the knowledge of 57 Fe density is taken 

39 into account. 

40 The Lagrangian that couples axions to nucleons is: 



C = aijji>y 5 (g (3 + g3T 3 )ip. 



(1.5) 
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Table 1. Detection Rates for TeC>2 (counts-kg 1 -d 1 ) 



S\/ a (GeV) 3.0xl0 5 5.0xl0 5 7.0xl0 5 9.0xl0 5 1.1x10' 



0.15 


0.070 


0.014 


0.007 


0.00273 


0.00122 


0.20 


0.349 


0.069 


0.037 


0.0136 


0.00609 


0.25 


0.840 


0.166 


0.089 


0.0327 


0.01466 


0.30 


1.542 


0.304 


0.164 


0.0601 


0.02694 


0.35 


2.457 


0.485 


0.262 


0.0958 


0.04293 


0.40 


3.585 


0.707 


0.382 


0.140 


0.06263 


0.45 


4.924 


0.972 


0.525 


0.192 


0.08604 


0.50 


6.476 


1.278 


0.690 


0.253 


0.11316 


0.55 


8.241 


1.627 


0.878 


0.321 


0.14398 



41 Here, go and 53 are the iso-scalar and iso-vector coupling constants, and T3 is a Pauli matrix. 

42 The axion-nucleon coupling constants are defined by [12]: 

o/6.2 x 10 6 GeV\ /3F-D + 2S\ , . 

,o = -7.8 x 10-( ) ( ) (1.6) 

o/6.2 x 10 6 GeV\ /,„ ^l-z\ , „ N 

53 = -"xl0 -»( s ){(D + F)—) (1.7) 

43 where F~0.48 and D=0.77 are invariant matrix elements of the axial current [13]. The Spin- 

44 Muon Collaboration lists a range for S of 0.15 < S < 0.50 at 95% C.L [14] while Altarelli et. 

45 al give the range 0.37 < S < 0.53 [15]. Both go and 53 contribute to the axion branching ratio. 

46 This is critical to the formulation of the energy loss within the sun due to axion emission. 

47 Combined with the flux rates listed in Ref. [8] the axion detection rates were determined and 

48 listed Table 1. 



49 2 Experiment 

50 In the analysis shown in the present work we will focus on the results from a specific R&D 

51 run dedicated to the CUORE experiment. 

52 CUORE will be an array of 988 tellurium dioxide crystals, each with size 5x5x5 cm 3 

53 and weight 750 grams, which will be operated as bolometers at a temperature of about 10 mK 

54 to search for neutrino-less double beta decay of 130 Te and other rare events. A description of 

55 the CUORE technique and the basic principles behind bolometers is given in [16], [17], [18], 

56 and [19]. 

57 The crystals to be used for CUORE are produced at the Shanghai Institute for Ceramics 

58 of the Chinese Academy of Science (SICCAS), and are shipped to the Laboratori Nazionali 

59 del Gran Sasso (LNGS) located in Assergi, Italy. The crystals are produced and shipped 

60 in batches to LNGS. Four crystals are taken at random from each batch to measure their 

61 radioactive contamination levels and evaluate their performance at low temperature. Each 

62 one of these runs is called a Cuore Crystal Validation Run, or CCVR [20]. This paper will 

63 analyze data from the second run, known as CCVR2. 

64 The four Te02 crystals have a total active mass of 3 kg and are mounted in a specially 

65 designed copper frame which is placed inside a dilution refrigerator at LNGS. The cryostat 
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f a [GeV] 

Figure 2. Expected rate in the axion region as a function of the f a axion constant for different values 
of the nuclear S parameter. The horizontal line indicates the upper limit obtained in this work. 

66 is maintained at a working temperature of approximately 8-10 mK throughout the runs 

67 duration. 

68 Attached to each CCVR crystal are two neutron transmutation doped germanium semi- 

69 conductors, NTDs. CCVR2 collected data for a total of 19.4 days for a total of 43.65 kg-d. 

70 Calibrations were performed in the middle and end of the run with a 232 Th 7 source inserted 

71 inside the Pb shielding, close to the cryostat outer vacuum vessel. 

72 Previous CUORE R&D experiments, such as CUORICINO [16] [21], had bolometers 

73 whose typical threshold was 50 keV. Since the Q- value for 130 Te is at 2527 keV [22] [23], a 

74 threshold at 50 keV was more than satisfactory for /3/3-decay searches. However, to investigate 

75 physical events at lower energies, new procedures were needed to lower the energy threshold. 

76 For this reason a special low-energy trigger was developed and applied offline, exploiting the 

77 fact that the standard CUORE DAQ used in the CCVR runs collects data continuously and 

78 with no hardware trigger using 125 Hz 18-bits digitizers [24]. 

79 This trigger algorithm maximizes the signal to noise ratio by filtering the data with a 
so known power spectrum and a reference signal shape. A full description of it is given in [25] . 
si More details can be found in [26]. 
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82 3 Results 



83 In order to reject thermal and microphonic noise, the pulses are selected by applying cuts on 

84 the pulse shape parameters. The efficiency of these cuts must therefore be taken into account 

85 in the analysis. As usually done in Cuoricino and in all CCVR runs, a calibrated current 

86 pulse is sent every 30 minutes to a small resistor (heater) that is glued on each crystal. This 

87 calibration pulse is used to correct for temperature variations in the cryostat that change the 

88 response of the bolometer and therefore its signal amplitude. It is also used to measure the 

89 detection efficiency of the crystals. 

90 The detection efficiencies (erj) were measured using the heaters [27] [28] [29] through a 

91 dedicated measurement performed at the end of the data-taking. The equivalence of the 

92 bolometers response to particle and heater pulses was demonstrated via a MonteCarlo simu- 

93 lation [30] that allows to infer eo for the bolometer without the heater. Since the detection 

94 capability may vary with time, the fluctuation of the peak rate with time (11%) of the low- 

95 est energy 7 line (a 4.7 keV line described below) can be used as an (over) estimate of the 

96 systematic error on the detection efficiency [26]. 

97 The estimated values of e D are 0.91 ± 0.10 for Bl and 0.83 ± 0.09 for B2 and B3, 

98 dominated by the systematics. The pulse shape cut efficiencies have been evaluated on the 

99 4.7 keV peak and are equal to 1. 

100 The low energy spectrum, below 50 keV, is shown in Fig. 3. For this work we choose 

101 to include only the channels whose thresholds are lower than 5 keV as to include the peak 

102 at 4.7 keV. The physical interpretation of this peak is still being analyzed, but the lack of a 

103 full understanding of its origin has no consequence on the present work, because we used its 

104 width only to determine the energy resolution in the 14.4 KeV region. 

105 Energy calibration in this very low energy region is not done with the 232 Th lines, which 

106 are normally used for higher energy studies. The energy region between 2.5 and 300 keV is 

107 calibrated using a set of metastable Te lines that result from cosmogenic activation. The 

108 crystals spend a few weeks above ground while they are being shipped by sea 1 between the 

109 production site in Shangai, China and arrival at the underground storage site at LNGS. The 
no main lines used in the calibration for present work are reported in Table 2. 



Energy [keV] 


Source 


Life-Time (days) 


30.4912 


Sb X-ray 




88.26 ± 0.08 


127mrp e 


109 ± 2 


105.50 ± 0.05 


129mrp e 


33.6 ± 0.1 


144.78 ± 0.03 


125m'pg 


57.40 ± 0.15 


247.5 ± 0.2 


123mr|-ig 


119.7 ± 0.1 


293.98 ± 0.04 


121mrp e 


154 ± 7 



Table 2. List of 7 lines from meta-stable Te isotopes used in this analysis for the energy calibration 
in the energy region between 2.5 and 300 keV. The lines are available from cosmogenic activation of 
Te during shipment, and have half-lives spanning from 33.6 days and 119.7 days. 



1 The shipment of two (out of four) crystals was actually made by airplane, which induces a slightly higher 
activation. 
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Energy [keV] 

Figure 3. Energy spectrum for the low energy region between threshold and 40 keV. The two peaks 
at 4.7 keV and 30.49 keV are those used to study the energy resolution at low energy. The axion 
region around 14.4 keV is magnified in the inset with the fit result shown as red curve. No peak is 
observed at 14.4 keV, the Ml transition energy of 57 Fe for solar axions. 

in As a further check of the calibration, we faced to a bolometer a 55 Fe source. The X- rays 

ii2 produced, with nominal energy between 5.888 and 6.490 keV, were shifted by only 48±16 

in eV. 

ii4 The energy calibration is obtained from a third degree polynomial fit. The complete low 

us energy spectrum is modeled with two exponentials, one for the region lower than 5 keV, one 

116 for the background above 7 keV, and also a Gaussian indicative of the peak around 4.7 keV: 

f 9 lobal(E) = ax ■ e~^ E + a 2 ■ e~^ E + -^e^b^ . (3.1) 

V27TC7 

117 Extracting the fit parameters obtained using Eq. 3.1 applied to Fig. 3 yields 04.7 = 0.31 ±0.04 
us keV. This value is then applied to the fit in the axion region from 11 to 18 keV. It is not 

119 as reasonable to extract limits based on the overall global fit since the fluctuation in the 

120 background is more pronounced in the spectrum shown in inset of Fig. 3, the region where 

121 the 14.4 keV axion interaction would be observed. 

122 To this spectrum a "modified" fit is applied to extract limits on the axion's coupling 

123 constant and mass. First, the background needs to be modeled. The background function 

124 appears to be flat, but to ensure validity first, second and third order polynomials are tested 

125 along with an exponential. For each, a Pearson's-^ 2 test is performed to determine the most 
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fa X ion{E) = axe** + - 1 =^Exp v - z J (3.2) 

Lai 



accurate background fit. We find that an exponential models the background most accurately, 
so this is used with a Gaussian function to represent the 14.4 keV axion peak. The variance for 
the Gaussian will be the previously found 04.7. The bolometers resolution is not dominated 
by signal fluctuations and we know that the energy dependence of the detector resolution is 
weak; no significant change is expected from 4.7 keV to 14.4 keV 2 . The signal plus background 
model which we have used in the fit is therefore: 

Jj-E , Q 2 jj, \{E — Eua) 2 ' 

\/27r0"4.7 - ^"4.7 

No excess in the axion region is observed. We set therefore a lower bound on the axion 
decay constant f a based on the background fluctuations. 

There are 116±12 counts in the energy interval between [-E14.4 — v2o"4, 7 ,-©14. 4 + v2o"4.7] 5 
corresponding to 2.65±0.27 c-kg -d . 

Assuming Poisson statistics and including the combined effect of a fluctuation of the 
background and of the 15% systematic error induced by cut efficiencies and the uncertainty 
in the solar axion rate, we assert that the axion detection rate is < 0.63 c-kg^-d -1 at a 95 
% C.L. Using S=0.55, this places a bound on the axion coupling constant of f a > 0.76 x 10 6 
GeV at 95% C.L. 

The 95% C.L. and the 1-a bounds as a function of the nuclear S parameter are shown 
in Fig. 4. 



143 4 Summary and Conclusions 

144 An experimental search was performed for axions from the solar core from the 14.4 keV 

145 Ml ground-state nuclear transition in 57 Fe. The detection technique employed a search for 

146 a peak in the energy spectrum at 14.4-keV when the axion is absorbed by an electron via 

147 the axio-electric effect. The cross section for this process is proportional to the photo-electric 

148 absorption cross section for photons. In this pilot experiment 43.65 kg-d of data were analyzed 

149 resulting in a lower bound on the Peccei-Quinn energy scale of f a > 0.76 x 10 6 GeV for a 

150 value for the flavor-singlet axial vector matrix element of S=0.55; bounds are presented in a 

151 graph for values 0.15 < S < 0.55. 

152 With the numbers quoted in the text, the limit on f a translates into a mass limit m a < 6 

153 eV, significantly more stringent than in recent results obtained with 57 Fe detectors [31] and 

154 by the Borexino experiment [32] [33]. 

155 The CUORE experiment will have about 740 kg Te02- With live-time of 5 live years 

156 the exposure will be 1.4 x 10 6 kg-d. With a similar background as the one reported here the 

157 expected sensitivity on f a will be roughly increased by an order of magnitude. 
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Figure 4. Exclusion plots for one and two sigma confidence levels for f a as a function of the nuclear 
S parameter values (see Eq. 1.7 for the meaning of the S parameter). 
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